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Abstract 

We present predictions for a variety of single-inclusive observables that stem from the produc- 
tion of charm and bottom quark pairs at the 7 TeV LHC. They are obtained within the FONLL 
semi-analytical framework, and with two "Monte Carlo + NLO" approaches, MC@NLO and 
POWHEG. Results are given for final states and acceptance cuts that are as close as possible to 
those used by experimental collaborations and, where feasible, are compared to LHC data. 
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1 Introduction 



Measurements of charm and bottom production are among the most interesting studies of QCD 
dynamics that have emerged from the analysis of the first LHC data at 7 TeV centre-of-mass energy. 
The large total cross sections, expected to be of the order of 5 mb (charm) and 250 /ib (bottom), 
have provided abundant data samples already with the ~ 50 pb _1 collected in 2010. A precise 
knowledge of charm and bottom quark cross sections and distributions is important in order to 
assess the accuracy of QCD calculations, but also as a test of theoretical predictions susceptible of 
being used to estimate backgrounds to new physics searches. This is all the more true since, in the 
past, early experimental measurements of bottom production at the Fermilab Tevatron seemed to be 
significantly larger than QCD calculations [THS]. It took a while (see |10pilj for a review) before these 
discrepancies could be resolved, through improvements in the accuracy of both the experimental 
measurements |12H14j and the theoretical predictions. In particular, from the theoretical side, the 
introduction of the "Fixed Order + Next-to-Leading Log" (FONLL) framework |15yi6| has shown 
that the discrepancies are largely reduced if a consistent use is made of the fragmentation function 
information coming from e + e~ data. From the experimental side, the extension of the b production 
measurements to very small transverse momenta |12j has demonstrated good agreement with the 
fixed order QCD calculation [17j in a region where theoretical uncertainties due to fragmentation 
effects have very little relevance. 

These successful comparisons notwithstanding, independent checks at a larger centre-of-mass energy 
and with different measurements are of great interest. Firstly, Iow-pt production at the LHC 
energies probes values of the momentum fraction x smaller than at the Tevatron, and challenges 
QCD in a dynamical region where potentially large higher-order corrections need to be resummed. 
This is particularly true of production at large rapidity, which pushes x of one of the two initial- 
state partons to values of x ^ 10 -4 , and which can be studied by forward detectors like LHCb 
and the ALICE muon spectrometer. Secondly, the greater beam energy and the high luminosities 
of the LHC can push the kinematic reach to much larger transverse momenta, exposing another 
interesting dynamical regime, where the resummation of logarithms of big ratios (such as Pj-/ m Q, 

E^/rriQ, Pt/Pt i e ^ c -) ma y become crucial. 

In parallel with the most recent comparisons with Tevatron data mentioned above, new theoretical 
tools have become available. The successful matching of next-to-leading order QCD calculations 
with parton shower Monte Carlos (PSMCs) has led to the MC@NLO [TH1IT9] and the POWHEG 
[201ET] implementations, with matching to HERWIG (22j[23] and to PYTHIA [24], which allow one 
to obtain predictions for fully exclusive observables while guaranteeing that inclusive quantities 
retain full NLO accuracy. Predictions can now therefore be evaluated within different frameworks 
(see also |25y26j for an independent analysis) and compared among themselves as well as with the 
data. 

The purpose of this paper is threefold. Firstly, we compare theoretical predictions obtained within 
the FONLL, MC@NLO and POWHEG frameworks for realistic observables (heavy mesons/hadrons, 
leptons from heavy hadrons, J/V> from 6-hadron decays). Secondly, we compare these predictions, 
calculated within the cuts employed by the experimental analyses, with the available data. Thirdly, 
we provide a detailed public record for predictions that have been transmitted to the experimental 
community over the past two years, and that have been used in the comparisons with data that are 
documented in several experimental papers. 
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2 Description of the theoretical frameworks 



2.1 FONLL 

In the FONLL framework [15] , one matches fixed next-to- leading order (NLO) QCD [271128] with 
all-order resummation to next-to-leading log (NLL) accuracy in the limit where the transverse 
momentum (pr) of a heavy quark is much larger than its mass (m) [29]. It allows one to calculate 
predictions for one-particle inclusive distributions of a heavy quark (or heavy hadron), while the 
degrees of freedom of the other particles in the event are integrated ovei0- The accuracy of the 
FONLL calculation can be denoted as being NLO+NLL, with the understanding that the logarithms 
resummed up to next-to-leading accuracy are of the form a™ log™ (py/m) and a™ log n_1 (pT/m). The 
perimeter of the meaning of the name "FONLL" is often stretched to denote also the determination, 
from e + e~ data [30], of the non-perturbative fragmentation parameters necessary for a successful 
description of physical differential distributions. 

FONLL has been used extensively to predict bottom |16]I17] and charm [31] production data at the 
Tevatron and at RHIC [32]. In all cases satisfactory agreement between theory and data was found. 
The framework and the parameters employed in those predictions have not been modified since 
(except for a more systematic determination of the non-perturbative fragmentation parameters 
performed in [30], which confirmed earlier results). 

A prediction for a single inclusive distribution (typically transverse momentum (pr), rapidity (y) 
or pseudorapidity (77) ) of a particle I is obtained within FONLL as a numerical convolution of a 
perturbative cross section daq ONLL with a non-perturbative fragmentation function Dq^ Hq and 
possibly a decay function describing, for instance, the hadron weak decay into a lepton: 

1 FONLL rl^FONLL „ n iVP „ weak /-, \ 

The integral of the fragmentation functions Dq^ Hq , for a given heavy-flavoured hadron Hq , will 
be called the fragmentation fraction, f(Q—> Hq). The parameters for Dq^ Hq are best determined 
from e + e~ collisions data, and a systematic survey is given in [30| . Weak decay spectra g^^i and 
the corresponding branching ratios are similarly extracted from experimental data. 

All the FONLL predictions published in this paper, as well as others corresponding to different cuts, 
input parameters or distributions, can be obtained from a publicly accessible web page [33] . 



2.1.1 Non-perturbative fragmentation 

For completeness, we summarise here the parameters that we have used in this paper. 

For bottom production, the functional form chosen for the parametrization of the non perturbative 
fragmentation function is a Kartvelishvili et al. distribution [34] : 

Df4 = (« + l)(« + 2)/(l-x). (2) 

We choose mj = 4.75 GeV as central value for the bottom quark pole mass, and the range mj, = 
4.5 — 5 GeV for estimating the associated uncertainty. In this work we adopt the values of the 

x We note that this implies that this approach does not allow the study of correlations between the heavy quark 
and antiquark. 
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fragmentation parameter a obtained in ref. [30J, namely a = 24.2 for nib = 4.75 GeV, a = 26.7 
for nib = 4.5 GeV and a = 22.2 for nib = 5 GeVl These parameters were fitted in |30| to the 
LEP data relevant to the production of a mixture of 6-hadrons |35|I36]. No data are available for 
the individual hadrons, for instance or B°. One is therefore forced to assume a similar non- 
perturbative fragmentation for all of them. This assumption is likely fulfilled to a large extent, with 
the possible partial exception of A& production, which however only contributes a minor fraction 
of the 6-hadron yield. We shall collectively denote the 6-hadrons by Hb, and assume that each b 
quark eventually produces a 6-hadron, i.e. the fragmentation fraction f(b — > Hb) is equal to one 
(and equivalently for the b). The function in eq. ([2]) is normalised to one: when used for a specific 
state the appropriate fragmentation fraction (e.g. f(b — > B + )) must additionally be provided. 

For charm production the situation is more complex. On one hand, experimental data are avail- 
able for individual mesons (D*, D^, D° and D°, D s ). On the other, at least partial theoretical 
understanding exists on the differences and similarities in the fragmentation of a heavy quark into 
a pseudoscalar (D) or a vector (D*) mesons (see e.g. [37]). This understanding can therefore be 
exploited in order to minimize the number of parameters to be extracted from experimental data. 
To this aim, in ref. |31j the non-perturbative fragmentation functions into different charmed mesons 
were constructed exclusively in terms of the c — > D* fragmentation, whose single parameter was 
extracted from ALEPH e + e~ data [38]. The decays of D* into D states was modeled theoretically, 
and the various branching ratios were extracted from data. Primary D production from c frag- 
mentation was described in terms of the same non-perturbative parameter fitted to c — > D* data, 
though a different functional form, as computed in [37J, was used. The detailed results for the 
non-perturbative fragmentation functions D^ D ,, D^ D+ and D^ D0 are given in eqs. (10), (9) 
and (5) of [H] respectively, and we do not repeat them here for brevity. They depend (besides the 
branching ratios of experimental origin) on a single non-perturbative parameter r determined from 
data. We used the values r = 0.1 with m c = 1.5 GeV, r = 0.06 with m c = 1.3 GeV and r = 0.135 
with m c = 1.7 GeV. These values also correspond to a A^ = 5 fit as detailed in [30] . 

2.1.2 Theoretical uncertainties 

The 'central' FONLL prediction is computed by setting the renormalisation and factorisation scales 
equal to the transverse mass, Hr : f = A*o — \Jv\ + m2 ' nere > rn and pt are the heavy quark mass 
and transverse momentum respectively. The theoretical uncertainty is estimated as a combination 
of factorisation and renormalisation scale variations, heavy quark mass variation, and uncertainty 
associated with Parton Distribution Functions (PDFs). The three uncertainties will normally be 
combined in quadrature. 

• In order to avoid accidental compensation between the fj,p an d the [Ir dependence of the cross 
section, which may occur if the two scales are set equal, we compute the scale uncertainty by 

2 These parameters differ from those employed in |16II17|[32] . where a = 29.1 was used as a central value with 
■nib = 4.75 GeV, together with a = 34 with mj = 4.5 GeV and a — 25.6 with mj = 5 GeV (note that in the text of 
refs. [171132] the values for a corresponding to rxib = 4.5 and mj, = 5 GeV are unfortunately inverted, although the 
numerical results were obtained with the correct ones). The choice made in 16, 17,32 corresponded to using a fit to 
the N — 2 Mellin moment of the e + e~ data, whereas the values used in this paper correspond to a fit to the N = 5 
moment. The reason for a different choice is that the N = 5 moment is more appropriate for describing steeply falling 
Pt distributions like the ones found at large pr at the LHC. One should note, however, that switching from the old 
parameters set to the new one can be seen to lower the predictions by 5-10% at most: it does not affect therefore the 
good agreement with the data found in the past. 
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varying them independently over the range 0.5 < £rf — 2, with the constraint 0.5 < £r/£f < 
2, where ^r^f = Mh,f/a*o- I n practice, the cross section is evaluated (using the central mass 
value and PDF set) in the seven points 

(6i,£f) E {(1,1), (0.5,0.5), (2,2), (0.5,1), (1,0.5), (2,1), (1,2)}, (3) 

and the envelope is taken. This envelope defines, at each point in the distribution one is 
considering, the two extremes of 

d(j +A + , scales ,^ 



• The three mass values m = 1.5, 1.3 and 1.7 GeV are used for charm, and m = 4.75, 4.5 and 5 
GeV for bottom. Non-perturbative parameters are adjusted for each mass to their appropriate 
valuejl and the cross section is evaluated setting = 1 and with the central PDF set. The 
envelope is then taken, defining 

d(J -A ■ ( 5 ) 

*-» — ,mass v ' 

• The PDF uncertainty, where available, is evaluated (setting £n t F = 1, and the mass equal to 
the central value) as suggested by the specific PDF set used (see e.g. [39]), yielding 

d --A™ • (6) 
Our default PDF set will be CTEQ6.6 [3D], unless otherwise stated. 

The full uncertainty band of the FONLL predictions will be given by 



da + _i + _ (7) 



with 



,PDF (8) 



2.2 NLO+PS approaches 

Heavy flavour production has been available for quite some time in NLO+PS (Next-to-leading order 
plus parton shower) models, namely in MC@NLO [19] and in POWHEG [21j . These programs, in 
conjunction with existing parton shower programs, are capable of generating fully exclusive final 
states, maintaining the next-to-leading order accuracy for inclusive observables. 

We remark that, unlike the FONLL approach, these generators resum only a subset of all con- 
tributions enhanced by large logarithms of the ratio of the transverse momentum of the heavy 
quark over its mass. In particular, the so called gluon splitting and flavour excitation production 
mechanisms are implemented only at order 0(a%), i.e. the lowest order contributions to these pro- 
cesses. In contrast, FONLL includes for these contributions all terms of order a| x (alogpT /Tn) n 
and a s x (as log pT/m) n . On the other hand, at small and moderate transverse momenta, the 
NLO+PS approaches are superior to FONLL, since they have the same accuracy in this region but, 
unlike FONLL, provide a complete and fully exclusive description of the final state. 

3 Note that this adjustment translates into a non negligible reduction of the sensitivity of the hadron-level cross 
section on the heavy quark mass value at large transverse momentum, as one may expect from the fact that neither 
the heavy quark mass nor the non-perturbative fragmentation are physical observables, and therefore their variations 
must compensate each other in their interplay. 
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In the present comparison, the NLO+PS methods should be viewed as approaches that work at 
small and moderate energies, and that are bound to fail at some large transverse momentum scale. 
To some extent, the present work helps to assess the range of validity of these approaches. 

It should be also kept in mind that the non-perturbative part of the shower that leads to the 
formation of the b hadrons is handled by the Parton Shower program alone. The corresponding 
parameters are tuned using final-state observables reconstructed with particles emerging from the 
parton shower, but with the hard production cross section of the shower Monte Carlo, i.e. not that 
of the next-to-leading order ones used in the NLO+PS methods. This is in contrast with what is 
done in FONLL, where use is made of a non-perturbative fragmentation function which is fitted to 
e + e~ data with a theoretical calculation based on the very same underlying FONLL approach. 

We point out that the parameters that control cluster decays in Fortran HERWIG [41J (PSPLT, 
CLDIR, and CLSMR) may be assigned (depending on the tuning adopted) specific values for ^-flavoured 
clusters. On the other hand, c-flavoured cluster decays are treated in the same manner as light quark 
ones. Although one can envisage to relax this constraint, and to introduce and tune parameters 
relevant to c-flavoured clusters only, such a possibility is not given in the official HERWIG versions, 
and has not been considered here. 

2.2.1 MC@NLO 

The MC@NLO formalism has been introduced in ref. [18], and aims at a consistent matching between 
NLO QCD corrections for a given process, and parton showers. The relevant technical details are 
given in the quoted references, and we shall omit them here. We limit ourselves to recall that 
in the context of MCQNLO the matching prescription amounts to modifying the short-distance 
cross sections relevant to the NLO computation, by including the so-called Monte Carlo (MC) 
subtraction terms, that are responsible for removing any double counting at the NLO. The MC 
subtraction terms can be computed in a process-independent manner, but they are still dependent 
on the particular PSMC one adopts for the shower phase. In other words, each PSMC requires 
a set of MC subtraction terms, which can be obtained by formally expanding the PSMC results 
to the same order in a s as the corresponding NLO contribution to the parton-level cross section 
(i.e. that of the real-emission matrix elements). Furthermore, their structures are such that all 
non trivial process-specific information is contained in the Born matrix elements. These matrix 
elements are multiplied by kernels whose analytic forms depend solely on the shower variables used 
by the PSMC to generate the elementary branchings, and on the identities of the partons involved 
in such branchings. 

Although ref. [18] formulated the solution of NLO+PS matching in general terms, practical applica- 
tions there and in subsequent papers have been restricted to the choice of Fortran HERWIG [41H43] 
as PSMC. Recently, MC subtraction terms have been computed which are relevant to the matching 
with Herwig++ [44] and, for processes that feature only initial-state emissions, with PYTHIA 6.4 [24] 
(see refs. ^ andP] respective* . In the present work, the MC@NLO results have been obtained 
with Fortran HERWIG (v6.520, which, for the processes considered here, gives results identical to 
those of v6.510), using the implementation of QQ production presented in ref. [19]. We point out 
that, although such an implementation technically only requires the mass of the heavy quark to 
be larger than Aqcdj the conservative choice has been made of not including cc production in the 

4 The matching with PYTHIA (including the p T -ordered shower versions) for a generic process is being achieved 
in the context of aMC@NLO - see http://amcatnlo.cern.ch for more information on this project. 
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public MC@NLO package. This is due to the fact that, because of the relative lightness of the charm 
quark, contributions may be important where a light-parton scattering such as gg — > gg initiates 
showers that eventually feature a g — > cc branching. These kind of contributions are not included in 
an NLO+PS matched computation, unless the aforementioned branching is the first occurring in a 
given shower. Although this mechanism is relevant also in the case of bb production, it is expected 
to be less important there than in the case of cc final states. It is therefore interesting to com- 
pare NLO+PS predictions with FONLL, which does take correctly into account the contributions 
discussed above. The default scale choice for QQ production in MC@NLO is 

= g ( m T(Q) + m 2 T {Q)) , ml=pl + m 2 , (9) 

where p T is the transverse momentum w.r.t. the beam line (note that, in the real-emission contri- 
butions to the NLO cross section, the transverse momentum of the Q is generally not equal to that 
of the Q). 

2.2.2 POWHEG 

The POWHEG implementation of heavy flavour production is described in detail in ref. |21j . It 
is now available at the POWHEG BOX website, http://powhegbox.mib.infn.it/. It can be used to 
generate events with either tt, bb or cc pairs, and it is based upon the heavy flavour production 
next-to-leading order calculation of refs. [2"7" l l28 l l4"T] . The output is an event file (in the Les Houches 
format |48| ) . that can be fed through any Shower Monte Carlo program that complies with the 
requirements of the Les Houches Interface for User Processes (LHIUP), like the Fortran and C++ 
versions of PYTHIA and HERWIG, in order to generate complete events. 

The default scale choice used in this implementation is given by \jp\ + m2 > where px is the trans- 
verse momentum of the heavy flavour in the underlying Born configuration. This is the structure 
of the QQ event before the hardest radiation has been generated, as defined in refs. [39l[50]. It is 
obtained as follows. One performs a longitudinal boost of the heavy quark pair system such that 
the system has zero rapidity after the boost, then a transverse boost is performed such that the 
transverse momentum of the system vanishes, and then the inverse of the initial longitudinal boost 
is performed. Thus, although the total cross section in POWHEG corresponds to the NLO result, 
as in MC@NLO, one cannot expect the two generators to yield identical cross sections because of 
the different scale choice. 

POWHEG has been interfaced to PYTHIA version 6.4.25, and HERWIG 6.510, keeping always the 
default values of the parameters in the Monte Carlo. A non negligible sensitivity to the Monte Carlo 
tune is expected for the observables considered in this work. An example is given in figure [TJ where 
D + and B + production as predicted by POWHEG with PYTHIA is studied for different PYTHIA 
tunes in \y\ < 0.5. At large rapidities the sensitivity to the PYTHIA tune is similar, if not larger. 
A complete study of the tune dependence is outside the scope of the present paper, and all Monte 
Carlo results presented in this paper will employ the default tune, unless explicitly stated. 
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Figure 1: Sensitivity of POWHEG with PYTHIA predictions to the tune used in PYTHIA for D + 
(left plot) and B + (right plot) transverse momentum distributions in the central rapidity region. 
At large rapidities, we observe effects of similar size or larger. 



3 Numerical predictions and comparisons 
3.1 Open charm production 

We begin by considering the production of D mesons. The FONLL framework and its non- 
perturbative parameters are described in Section 12.1.11 In the following we concentrate on the 
production of D + mesons (the others being largely similar) and we use a value of the fragmentation 
fraction of /(c -)• D + ) = 0.238 with FONLL. The left plot of Fig. [2]shows the central FONLL predic- 
tion, compared with MC@NLO, POWHEG with PYTHIA (denoted by POWHEG-PY) and POWHEG 
with HERWIG (denoted by POWHEG-HW). The right plot shows the scale and PDF contributions 
to the FONLL uncertainty, as well as the overall systematics. The ALICE data from ref. |51j are 
also compared to the theoretical prediction, showing a good agreement within the large systematics. 
We note that all four theoretical predictions are in fair agreement with each other within the overall 
uncertainty. The large dependence, in shape as well as overall normalisation, of the Monte Carlo re- 
sults on the specific tune, as shown in section [2.2.2l should also be kept in mind. Differences appear, 
on the other hand, in the forward region 4.5 < \y\ < 5, as shown in Fig. [31 POWHEG-PY agrees 
well with FONLL, while the shower evolution carried out with HERWIG (whether in the MCQNLO 
or in the POWHEG framework) leads to a harder spectrum. These differences largely exceed the 
systematics quoted from the FONLL calculation, in particular if we consider that the sightly bigger 
uncertainty at large px is mostly due to the PDFs, and is therefore entirely correlated among the 
different predictions. It will be very interesting to see the first data on forward D + production at 
large pt from LHCb. 

Additional data on D° and D* + production in the central region, and comparisons with the FONLL 
predictions, are reported by the ALICE Collaboration in Fig. 5 of |51| . Similar data, at y/S = 
2.76 TeV, are reported in [52]. Preliminary data from ATLAS are also available, and compared to 
FONLL in Fig. 1 and 2 of [53J. In either case a fair agreement is found, with the data mostly centred 
on the upper edge of the theoretical uncertainty band, mirroring quite closely the comparisons with 
the Tevatron D mesons data [31 j . This may indicate a preference of the data for a value of the 
charm mass smaller than our default of 1.5 GeV. 
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Figure 2: Transverse momentum distribution of D + mesons at central rapidity, \y\ < 0.5. Left plot: 
comparison among the central predictions of our four benchmark calculations, FONLL, MC@NLO, 
and POWHEG with PYTHIA or HERWIG showers. Right plot: theoretical systematics for the 
FONLL calculation, and the comparison with data from ALICE [51]. For the systematics we show 
the individual scale and PDF components, as well as the combined total (which includes mass 
variation, as described in the text). 



3.2 Open bottom production from inclusive and fully-reconstructed H b — ^ D + X 
decays 

The theoretical predictions for central production of B + mesons {\y\ < 0.5) are presented in Fig. [H 
The quality of the agreement among the various predictions is similar to the one seen in the D + 
case above. Other features, like - at large rapidity - the progressively larger theoretical uncertainty 
as a consequence of a larger PDF uncertainty, or the larger variance between the NLO+PS and 
the FONLL predictions (see Fig. [5]), are also similar. The potential sensitivity of the Monte Carlo 
results to the specific tune used (though smaller than in the D-mesons case, see figure [1]), should 
be kept in mind in this case too. 

The first measurement of 6-hadron production at the LHC was performed by the LHCb Collab- 
oration |56j . The pseudorapidity distribution in the region 2 < r/ < 6 was shown to be in good 
agreement with NLO and FONLL predictions (see Fig. 5 of [56JL The measured total cross section 
in this region was found to be (averaging over b and b hadronsjj 

a LRCh (pp -> H b , 2 < 7] < 6) = 75.3 ± 5.4 ± 10.0 (ib (10) 

5 Note that, as explained in |56| . the measurement can change to 89.6 ± 6.4 ± 15.5 fib if 6-hadron fractions measured 
at the Tevatron rather than those measured at LEP are used in converting the number of events to an Hi, cross 
section. 
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to be compared with the FONLL predictiorjl 



a FONLL (pp -> H b , 2 < 7] < 6) = 70.8 tfd ^ (11) 

As remarked above, the contribution of 6-baryons are included in this theoretical value assuming 
they fragment like the B mesons. 

More recently, LHCb has published the measurement of the production cross section in the 
2 < y < 4.5 rapidity region [57], using fully reconstructed B^ — > J/ipK^ 1 decays: 

a Lnch (pp -> B ± , < p T < 40 GeV, 2 < y < 4.5) = 41.4 ± 1.5 ± 3.1 fib (12) 

to be compared with the FONLL prediction (which includes a fragmentation fraction f(b — > B~) = 
0.403) 

a FONLL (pp -> 5 ± , < p T < 40 GeV, 2 < y < 4.5) = 40.1 +H° S fib . (13) 

Good agreement in seen also in the differential p B± spectra, in the range up to 40 GeV, as shown 
in Fig. 2 of [57]. 

The CMS collaboration has published measurements of transverse momentum and rapidity distri- 
butions of B + [53], B° [55] and B® [58] mesons. They report, for the total visible cross sections at 
7 TeV LHC, 

a CMS (pp -> B°, > 5 GeV, \y B \ < 2.2) = 33.2 ± 2.5 ± 3.5 /x& (14) 
a CMS (pp -> £ + , pf > 5 GeV, < 2.4) = 28.1 ± 2.4 ± 2.0 ± 3.1 //& (15) 
cj CMS (pp 5°, 8 < < 50 GeV, < 2.4) x BR(5 S ° -> J/00) = 6.9 ± 0.6 ± 0.6 raft, (16) 



6 Here and elsewhere, where the quoted FONLL predictions already appeared in experimental papers. They were 
originally calculated according to the framework defined in this paper, and provided as private communications to 
the experimental collaborations. 
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Figure 4: Transverse momentum distribution of B + mesons at central rapidity, |y| < 0.5. Left plot: 
comparison among the central predictions of our four benchmark calculations, FONLL, MCQNLO, 
and POWHEG with PYLHIA or HERWIG showers. Right plot: theoretical systematics for the FONLL 
calculation, and the comparison with data from CMS [54|]55], rescaled to the \y B \ < 0.5 region. For 
the systematics we show the individual scale and PDF components, as well as the combined total. 



to be compared with the FONLL predictions 



a FONLL (pp -»• B°, p| > 5 GeV, \y B \ < 2.2) = 25.5 l^ 5 /j,b (17) 
a FONLL {pp -> B+, p$ > 5 GeV, \y B \ < 2.4) = 27.2 ^ 5 2 /jib (18) 
a FONLL {pp -> B°, 8 < Vt < 50 GeV - \V B \ < 2A ) x BR(B° -»■ J/ip<f>) = 4.5 ±tl ± !- 6 n6.(19) 

The FONLL uncertainties in eqs. (|17p - (|19p are due to renormalisation and factorisation scales, 
heavy quark masses and PDF, as detailed in section l2.1.2l Of these three sources, the first is largely 
dominant. In the FONLL predictions the fragmentation fractions f(b — > B°) = f(b — > B~) = 0.403 
and f(b -»■ B°) = 0.11 [59] have been included. These values, leading to a ratio f s /fd = 0.273, 
are consistent with the recent direct measurement of this ratio by LHCb [60]: f s / fd = 0-267^q'q2o- 
Additionally, BR(£>^ — > J/tp (ft) = (1.4 ±0.5) x 10 -3 [59] has been used, and the second uncertainty 
in the FONLL prediction of eq. (fl~9|) for B® production reflects the large uncertainty of this measured 
branching ratio. 

A good agreement with the experimental measurements within the respective uncertainties can be 
observed, confirming the latest comparisons between theory and Tevatron data which showed no 
significant excess in bottom hadroproduction compared to theoretical predictions. The pt spectra 
measured by CMS are compared to various predictions in Fig. 2 of ref. [55] and Fig. 2 of ref. |54j . 
The comparison of these data, rescaled to the rapidity region |y| < 0.5, with the FONLL predictions 
is shown in the right panel of Fig. |U 
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Figure 5: Same theoretical distributions as in Fig. HI for the forward rapidity region 4.5 < |y| < 5. 



3.3 Open bottom production from inclusive Hf, — > J/tp and — > ip(2S) decays 

The predictions for non-prompt J/tp and ip(2S), i.e. coming from 6-hadron decays, is obtained in 
FONLL by convoluting the distribution for 6-hadron production with a phenomenological spectrum, 
obtained from experimental data [HI] (that describes the momentum distribution of the J/tp and 
the tp{2S) in i?-meson decays), and by multiplying it by the appropriate branching fraction (again 
obtained from experimental measurements). In using decay spectra from B meson decays and 
applying them to all 6-hadron decays we make the reasonable assumption that a potential difference 
in the small fraction (~ 10-20%) of non- meson 6-hadrons is going to be inconsequential. Moreover, 
it should be noted that strictly speaking ref. [61] gives the momentum of the quarkonium in the 
T(4(S I ) rest frame rather than in the 5-meson one. The difference has a small rms spread (0.12 
GeV according to [61]) and is not expected to affect significantly the final quarkonium spectrum 
in our convolution. In the following we shall use the branching ratio BR(6 — > h), defined as 
BR(6 — > h) = Y2h /(^ ~~ H-b) x BR(-ff& — > h + X), summed over all relevant b-hadrons H^. 

The fit to the B — > J/if; decay spectrum from Fig. 6 of [61] had been performed in preparation 
of ref. [17j . where non-prompt J/ip distribution in pp collisions at the Fermilab Tevatron had been 
calculated, and found in good agreement with the data, as well as with the MC@NLO predictions 
(see Fig. 5 of [32]). This paper uses the same spectrum. The fit to the tp(2S) spectrum from Fig. 8 
of |61] has instead been performed for this paper, and predictions for the non-prompt ip(2S) within 
FONLL are given here for the first time. 

Figure [6] compares predictions for the transverse momentum distributions for the non-prompt J ftp 
and tp(2S) production in the central rapidity region \y\ < 2.4. In this case the prediction from 
the POWHEG-PY implementation differ quite markedly from the others, and data could provide 
discriminating power. 

The first experimental measurements at 7 TeV from CMS [62,63j, LHCb [63] and ATLAS [65J have 
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Figure 6: Transverse momentum distribution of non-prompt J/ip in the central rapidity region 
\y\ < 2.4. Left plot: comparison among the central predictions of our four benchmark calculations, 
FONLL, MC@NLO, POWHEG-HW and POWHEG-PY. Right plot: theoretical systematics for the 
FONLL calculation. 



been compared to FONLL predictions and generally found in fairly good agreement. For example, 
LHCb |64j . after separating the prompt and 6-decay contributions, measures: 

cr LHCb (J/0 from H b ,p T < 14 GeV, 2 < y < 4.5) = 1.14 ± 0.01 ± 0.16 fib (20) 

to be compared with the FONLL prediction (which includes a branching fraction BR(6 —> J/ip = 
0.0116) 

a FONLL (J/0 from H b ,p T < 14 GeV, 2 < y < 4.5) = 1.16 t° H fib (21) 

The FONLL predictions describe also very well the J/0 pt spectrum, measured in the range 0- 
13 GeV, as shown in Fig. 9 of [53] . 

More recently, ALICE [66] reported: 

ct ALICE (J/0 from H b:PT > 1.3 GeV, \y\ < 0.9) = 1.26 ± 0.33 ±%f 8 fib (22) 
to be compared with the FONLL prediction 

a FONLL (J/0 from H b ,p T > 1.3 GeV, \y\ < 0.9) = 1.33 l{J;gj fib (23) 

A possible exception to the generally good agreement is the observation, made by CMS in [63J, that 
the experimental cross section for non-prompt J/ip and ijj{2S) production tends to fall off at large 
transverse momentum slightly faster than the FONLL prediction. More data at even larger px will 
help clarify this issue. 

In the same paper, CMS note that the measured non-prompt ip(2S) production is quite smaller 
than the FONLL prediction evaluated using the central value of the branching fraction BR(6 — >• 
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Figure 7: FONLL prediction of the ratio of the cross section for non-prompt (i.e. from 6-hadrons) 
production of ip{2S) and J/ip (each multiplied by its own branching ratio into muons) as a function 
of their transverse momentum, compared to experimental data from CDF [12] at the Tevatron and 
CMS [63] and LHCb [67J at the LHC. 
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ip(2S)) = 4.8 x 10 3 [59J. In order to better understand this potential discrepancy, it is useful to 
consider the ratio 



Rms), j/tp) = dpT i - r """JZ,7r. y (24) 



as a function of the transverse momentum of the quarkonium. Figure [7] shows the FONLL predictions 
for this ratio, compared to data from CDF at the Tevatron [12] and CMS [63] and LHCb [67] at 
the LHC. One can see that the agreement with the Tevatron data is acceptable when one uses 
BR(5 -> ip(2S)) = (4.8 ± 2.4) x 10^ 3 [59] and considers its large (O(50%)) uncertainty, which is 
not shown in the plot. The data from the CMS and the LHCb collaborations seem instead to be 
somewhat at variance with the CDF ones whereas, according to the FONLL prediction, no large 
difference should be present. One could still exploit the large uncertainty of BR(6 — > ip(2S)) to 
lower significantly the theoretical prediction and make it compatible with all measurements within 
the uncertainties. In fact, CMS have used their measurement of this ratio and the comparison to 
theoretical predictions to extract [63] a new value for the BR(6 — > ■0(2/5)) branching fraction, 

BR(6 -> ij(2S)) CMS = (3.08 ± 0.12 ± 0.13 ± 0.42) x 10~ 3 . (25) 

LHCb have more recently performed a similar extraction of BR(6 — > ip(2S)) from their data, and 
have obtained 

BR(6 -> 0(2S')) LHCb = (2.73 ± 0.06 ± 0.16 ± 0.24) x 10" 3 , (26) 
fully compatible with the CMS determination. 

When these new branching ratios are used in the FONLL prediction, instead of the (4.8 ±2.4) x 10~ 3 
quoted above, the results are, essentially by construction of course, in full agreement with the 
LHC measurements, as shown in Fig. UJ On the other hand, the older CDF data now appear to 
stand out, albeit being potentially marginally compatible with the LHC ones within the respective 
experimental uncertainties. Since the theoretical predictions for this ratio are largely insensitive 
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Figure 8: Transverse momentum distribution of leptons from heavy hadrons in the central rapidity 
region \y\ < 2.5. The curves represent the sum of the three processes H c — > I, Hf, — )■ £ and 
Hb — > H c — > i. No PDF uncertainties are included in the NLO band. Data from [TO] are also 
shown. Note that both muon charge states are included. 



to the specific experimental setup, be it the center of mass energy (Tevatron or LHC) or even the 
rapidity acceptance cuts, the discrepancy is likely of experimental origin. 



4 Open charm and bottom production from inclusive semileptonic 
decays 

In a fashion similar to that employed to describe J/ip and ip(2S) production from 6-hadron decays, 
one can describe the production of a lepton (e or fi) originating from the electroweak decay of 
a charm or a bottom hadron. As before, the decay spectrum and the branching fraction can be 
extracted from experimental measurements. This approach has been used with FONLL for the first 
time in |32j . where the prediction for the production rate of electrons from charm and bottom in 
pp collisions at RHIC was calculated. The results were eventually found in good agreement with 
measurements by the PHENIX |68| and the STAR [69J collaborations. 

Three separate processes contribute to the final yield of leptons i from a heavy hadron: H c — > £, 
Hfj — > I, and the secondary decay — > H c — > t. The first contributes mainly at small transverse 
momentum (pj- < 5-10 GeV), the second dominates at larger pt, while the secondary decay is 
largely negligible (see e.g. Fig. 3 of [321). All three processes have been modeled for [32], and we 
use here exactly the same setup. 

Figure [8] shows the cross section predicted by FONLL for the sum of the three processes, as they are 
largely indistinguishable experimentally. Two NLO-level predictions are also shown. Both make 
use of the fixed-order heavy quark production NLO calculation for hadronic collisions. In one case 
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Figure 9: CMS data for muon production from 6-hadrons decays compared to FONLL and POWHEG 
predictions. Note that both charge states (i.e. muons from both b and b) are included. 



(the upper band, blue in colour display, and labeled simply 'NLO') the same non-perturbative 
fragmentation functions also employed in the FONLL case are used. In the case of the magenta 
dashed curve instead, labeled 'NLO with NLO NP fit', the non-perturbative fragmentation func- 
tions have been fitted to e + e" LEP data for B and D* production using a fixed next-to-leading 
order calculation. The FONLL and the 'NLO' predictions differ significantly at large pt when the 
same non-perturbative fragmentation functions are used, a consequence of the resummation effects 
included in FONLL . This difference is however largely compensated - as expected - by the adjust- 
ment of the non-perturbative contribution, and one can see that the magenta dashed curve 'NLO 
with NLO NP fit' is indeed much closer to the FONLL prediction. It should be noted, however, 
that this prediction employing the NLO fit to LEP data is expected to be valid only in a limited 
region, where transverse momenta are comparable to the scale \x at which the fit has been performed 
(/i = M z ~90 GeV). 

Data for muon production are available from ATLAS [70] up to very large transverse momentum, 
and one can see in Fig. [8] how they seem to be better described by the resummed FONLL prediction 
than by the fixed order one. Comparing to 'NLO with NLO NP fit' rather than to 'NLO' may 
appear to largely wash out the edge of FONLL at large pr, but the agreement with the data seems 
to deteriorate at moderate pt, as it may be expected from the use of non-perturbative fragmentation 
functions fitted with a fixed order calculation and at the much larger scale given by the Z mass. 
In fact, the good agreement with the 'NLO with NLO NP fit' is obtained by trading away the 
correct description of multiple quasi-colliner emissions via resummation in exchange of their effective 
inclusion in the non-perturbative fragmentation functions fitted to the NLO calculation. While 
theoretically less ambitious, the latter procedure is also bound to fail when scales much different 
than Mz are probed. 

The ALICE Collaboration has also produced comparisons of distributions of leptons from heavy 
hadrons with FONLL predictions, i.e. muons in the forward rapidity region 2.5 < y < 4 (see for 
instance Fig. 3 of [71]) and electrons in the central region \y\ < 0.5 (see Fig. 11 of |72j): the overall 
agreement is similar to the one with the ATLAS data shown in Fig. [8j 

CMS also published results for transverse momentum and rapidity distributions for muons from 
6-hadrons |73| . The plots in Fig. [9] compare their results with FONLL and POWHEG predictions, 
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and are found to be compatible within uncertainties. From their measurements CMS extract a total 
visible cross section for muons from 6-hadron decays, 

a CMS (pp -> H b + X -»■ // + X', p^>6 GeV, < 2.1) = 1.32 ± 0.01 ± 0.3 ± 0.15 fib (27) 

to be compared with the FONLL prediction (which uses the branching ratios BR(6 — > £) = 0.0108 
and BR(6 -»• c -»• £) = 0.096 [59]) 

a FONLL (pp -»• fli + X -> /x + X', p£ > 6 GeV, \y^\ < 2.1) = 0.855 j#fg fib (28) 

Note that both charge states are included in these results. 



Expt 


Observable (pt in GeV) 


cj cxp 


0.FONLL 


Comments 


1 


LHCb [55] 


<j{H b ,2< 77 < 6) 


75.3 ± 11.4 fib 


70.8 1^3.3 ^ b 


average 6 + 6 


2 


LHCb [57] 


a(B ± ,p T < 40,2 < y < 4.5) 


41.4 ± 3.4 /ib 


40.1 lit" M b 


B~) = 0.403 


3 


CMS 55, 


a(i? , pf > 5, < 2.2) 


33.2 ±4.3 fib 


25.5 t^ 5 /ib 


/(6 -> B°) = 0.403 


4 


CMS [54] 


(i(B+, pf > 5, < 2.4) 


28.1 ±4.4 fib 


27.2 t^ 5 2 /ib 


/(M S-) = 0.403 


5 


CMS [58] 


(t(B s °, 8 <pf < 50, < 2.4) 
xBR(B° -> J/V»c/») 


6.9 ±0.8 nb 


4.5 j£g nb 
(includes BR 
uncertainty) 


/(6 -> B°) = 0.11 
BR(B° -> J/ty> 0) = 
(1.4 ±0.5) x 10" 3 


6 


LHCb [64] 


<7(i/ b -> Jyty.p* < 14, 2 < ^ < 4.5) 


1.14 ± 0.16 fib 


1.16 t°;f 2 /ib 


BR(6 ->■ J/ty>) = 0.0116 


7 


ALICE 66, 


ct(# 6 -> Jyty.p* > 1.3,1^1 < 0.9) 


1.26 ±0.16 fib 


1.33 ±°;«g M b 


BR(6 ->■ J/V>) = 0.0116 


8 


CMS [73] 


cr(ff b -> /J, p£ > 6, |^| < 2.1) 


1.32 ±0.34 fib 


0.855 Mb 


BR(6 ->■ £) = 0.0108 
BR(6 -> c -> =0.096 



Table 1: Summary of various cross-section measurements, compared against the FONLL predictions. 
The numbers labeling each measurement refer to the entries in Fig. [101 



5 Conclusions 

We have presented theoretical predictions for heavy quark cross sections for various observable 
final states within realistic acceptance cuts. In a number of cases, we compared the predictions 
of FONLL, MC@NLO and POWHEG with HERWIG and/or PYTHIA. These predictions are also 
compared directly with experimental measurements at the LHC, and this paper is meant to provide 
a detailed description and a record of these theoretical results. 

The curves presented in this article, as well as others directly delivered to the experimental col- 
laborations, were mostly obtained before any data were available. They have not in the least been 
influenced by the LHC measurements, and are a direct extension to LHC energy of a framework 
previously built and validated using LEP and Tevatron data. 

One can note that, while FONLL, MC@NLO and POWHEG-based predictions generally agree with 
each other in the moderate transverse momentum regions, some differences can be observed at large 
Pt- The small-pT agreement is expected, as by construction all three frameworks must return the 
NLO cross section when integrated over the whole phase space, and this result will be dominated 
by the small-pT region. At large pr, instead, differences can arise. Part of them may be due to the 
Monte Carlo based approaches only resumming the quasi-collinear logarithms to 'almost' leading 
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Figure 10: Graphical display of the cross section results listed in Table [TJ We show the ratios 
of experimental data over the FONLL calculations, with the sizes of experimental and theoretical 
systematics. The numbers labeling each measurement refer to the entries in Table HJ 



log accuracy rather than next-to-leading one, but this is expected to be a minor effect. More im- 
portantly, non-perturbative fragmentation contributions that had been adjusted to data in PYTHIA 
and HERWIG are likely not fully appropriate to be used in a matched context where hard gluon 
radiation to next-to-leading order level is also included. These non-perturbative parametrisations 
should therefore be re-adjusted to e + e~ data within the same MC@NLO and POWHEG frameworks 
that are then going be used to produce predictions for hadronic collisions, analogously to what is 
done with FONLL. Lack of an appropriate tuning of the non-perturbative hadronization of heavy 
quarks could provide an explanation for the discrepancy between the fragmentation function of D* 
mesons inside jets, measured in Ref. [73], and the theoretical predictions, based on NLO PSMCs, 
shown in that paper. 

The agreement between the LHC experimental data we discussed and theoretical predictions is 
otherwise generally very good. We summarize the various cross sections reported in this paper 
in Table [J and in Fig. [TU1 The comparison of theory and data at small pt, as well as at large 
rapidity, does not point to the existence of a new dynamical regime for the production of heavy 
quarks (e.g. one dominated by small- £ effects), beyond what can be accounted for by fixed-order 
NLO calculations. The measurements at large transverse momentum, furthermore, seem to favour 
a description based on FONLL, underscoring the importance of resumming large logarithms of 

Prl m Q- 

Other sets of measurements not discussed in this paper, such as the Er spectrum and angular 
correlations of jets containing b quarks [75J7B], are consistent with these findings. Some discrepancies 
with MC@NLO, on the other hand, have been reported recently in a study of angular correlations 
between 6-hadron pairs |77| . 
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It is important to consider that, while the overall theoretical systematics is typically large, most of 
it is highly correlated in different measurements and in different kinematic regions. This is certainly 
the case of mass, fragmentation and PDF systematics, but it is also true of the scale systematics. 
When these correlations are taken into account, the agreement with the experimental data is even 
more remarkable, particularly considering the efforts that were required, during the first 10 years 
of measurements at the Tevatron, to reach such a level of consistency. 
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